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Abstract 

We investigate the effect of mechanical strain on the dynamics of thin M 0 S 2 nanodrum resonators. 
Using a piezoelectric crystal, compressive and tensile biaxial strain is induced in initially flat and 
buckled devices. In the flat device, we observe a remarkable strain-dependence of the resonance 
line width, while the change in the resonance frequency is relatively small. In the buckled device, 
the strain-dependence of the damping is less pronounced, and a clear hysteresis is observed. The 
experiment suggests that geometric imperfections, such as microscopic wrinkles, could play a role 
in the strong dissipation observed in nanoresonators fabricated from 2-D materials. 
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Nanomechanical resonators fabricated from 2-dimensional layered materials, such as 
graphene and M 0 S 2 , are known to exhibit low quality (Q-) factors at room temperature [T]- 
[Tj. The spectral Q-factor of these devices is orders of magnitude below the values that 
can be achieved with top-down fabricated devices, such as silicon nitride nanostrings [8]. 
Time-domain measurements on M 0 S 2 resonators with a thickness down to a single layer, 
revealed that the low spectral Q-factor is in agreement with the energy relaxation rates [7], 
indicating that the line-width is limited by dissipative processes. Although several 
mechanisms have been proposed for the high dissipation, such as clamping losses, surface 
effects, and energy leakage to other vibrational modes [U lUI-lTB] , the dominant mechanism 
responsible for the excessive dissipation is not identihed. 

It is well known that the Q-factor of top-down fabricated micro-and nanomechan¬ 
ical resonators (MEMS and NEMS resonators) can be increased by introducing tensile 
strain [HHIIl. It is explained by considering a complex elastic modulus, E = Ei -|- iE 2 , 
where the real part, Ei, corresponds to the Hooke's law spring constant, and the imaginary 
part, E 2 , gives rise to dissipation (energy loss). The intrinsic Q-factor can then be written 
as Q = E 1 /E 2 . Applying tension increases the real (conservative) part of the elasticity. 
This results in an increase in the resonance frequency, which is proportional to a/Ei, and an 
increased Q-factor, which is proportional to the resonance frequency by Q = fo/linewidth. 
Previous studies have shown that in MEMS devices the imaginary part of the elasticity 
can be assumed constant, i.e., independent of strain dSl UHl DS], and that the tensile 
strain enhances the Q-factor via the real part of the elasticity, leaving the line-width of 
the resonance peak virtually unaffected [la Eo]. Since strain engineering is commonly 
applied to realize MEMS resonators with high Q-factors, it is interesting to investigate the 
strain-dependence of the Q-factor of mechanical resonators from 2-D materials, which are a 
10 to 1000 times thinner, and exhibit very low Q-factors at room temperature. 

Here we study the strain-dependence of the resonant properties of M 0 S 2 nanodrum 
resonators. In contrast to tuning the strain by attracting the drum towards an electrostatic 
gate electrode m m, by applying a pressure difference [22], or by chemical modihca- 
tion [231 ED, we use a piezoelectric bender to introduce strain. This enables precise control 
over the strain, and allows one to study the drum dynamics without exerting out-of-plane 
forces that could affect the shape of the drum. Both tensile and compressive strain can be 
introduced. Two devices are considered: one that is initially flat, and one that is initially 
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buckled. In the flat device we observe a weak dependence of the resonance frequency on 
the strain, but a surprisingly strong strain-dependence of the line-width. This indicates 
that the tensile strain enhances the Q-factor via a reduction of the dissipative part of the 
elasticity, E 2 . This is in sharp contrast to top-down fabricated MEMS resonators in which 
the Q-factor enhances through an increase of real part of the elasticity, Ei. In the buckled 
device, the changes in the Q-factor are less pronounced. Here we observe hysteresis that 
could indicate a conformational change of the material, and possibly hints at the underlying 
process that causes the strain-dependent damping. 

To fabricate suspended M 0 S 2 resonators, we start with a 100 pm thin silicon wafer 
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FIG. 1: (a) Optical photograph of the device. The circles marked A and B indicate the considered 
devices. The diameter of the holes is 5 pm, and the thickness of the flake is 15 nm, which corresponds 
to 30 layers, (b) Setup for strain-tuning of mechanical resonators. A thin oxidized Si wafer 
containing the M 0 S 2 flakes is glued onto a sheet of piezoelectric material. PZT: Lead Zirconate 
Titanate piezoelectric material, (c) Interferometric displacement detector. HeNe: Helium-Neon 
laser; PD: photodiode; LD: laser diode, (d) Frequency response showing the lowest vibrational 
modes of drum A with Vp = 0. The relation between the resonance frequencies f 2 /fi ~ 2, which 
indicates the drum behaves mechanically like a plate. 


with a 285 nm thick layer of thermally grown silicon oxide. Thin Si wafers have a low 
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bending rigidity, and this enables the generation of signihcant mechanical strain. Circular 
holes are etched in the silicon oxide by conventional electron beam lithography and dry 
etching. M 0 S 2 flakes are mechanically exfoliated and deposited onto the substrate using 
a dry transfer method 12 a. Figure 1(a) shows the fabricated device; the diameter of the 
considered drums, marked A and B, is 5 /xm, and the thickness is 15 nm, which corresponds 
to ~ 30 layers. 

The wafer containing the drum resonators is hxed onto a commercially available Lead 
Zirconate Titanate (PZT) piezoelectric sheet, with electrodes on top and bottom. The 
wafer and the piezoelectric sheet form a bimorph structure, as is shown in Fig.1(b), which 
bends when an electric held, Vp, is applied across the piezoelectric sheet. Depending on the 
polarity, compressive or tensile strain is generated in the M 0 S 2 drum. The motion of the 
resonator is detected using an optical interferometer, shown schematically in Fig. 1(c). The 
suspended part of the M 0 S 2 hake forms the moving mirror, while the silicon substrate acts 
as the reference mirror. When probing the drum with a Helium-Neon laser, the intensity of 
the rehected optical signal is modulated by the position of the membrane P [SHZl [ 2 S], and 
detected using a photo-diode. To measure the frequency response of the drum, it is driven 
photo-thermally P using a laser diode (A = 405 nm) with an rf-modulated intensity. The 
measurements are performed at room temperature and at a pressure of 10“'^ mbar. 

Figure 1(d) shows a measured frequency response of drum A, when the voltage applied 
across the piezo is zero. The fundamental resonance frequency is detected at fg = 9.9 MHz, 
with a Q-factor Qq = 19. The second mode is observed at fi = 20.7 MHz, with a Q-factor 
Qi = 43. These low Q-factors are typical for mechanical resonators from 2-D materials at 
room temperature in vacuum PHElp. The ratio between the resonance frequencies for 
the lowest two vibrational modes, fi/fo ~ 2, is in agreement with a plate-like resonator. 
This is as expected, since the thickness of the resonator is beyond the membrane to plate 
cross-over, which occurs for M 0 S 2 at approximately hve layers [B]. In a plate-like resonator 
the restoring force arises mainly from the bending rigidity, whereas in a membrane it arises 
from the tension. 

We now measure the response of drum A while introducing strain by applying 
a voltage to the piezo. For each voltage the resonance frequency, the Q-factor, and the 
line-width are obtained from a harmonic oscillator ht. Figure 2 shows the result: the left 
column represents a compression cycle with Vp < 0, and the panel on the right a tensile 
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FIG. 2: Resonance frequency (a) and Q-factor (b) as a function of the voltage applied to the 
piezo. The panels on the left represent one cycle with compressive strain (Vp < 0), the panels on 
the right one cycle with tensile strain (Vp > 0). Multiple cycles were performed, and the (linear) 
trends of the resulting strain-dependence of the resonance frequency and Q-factor were calculated 
by least-squares fits to the data, and plotted in the histograms in the insets. Panel (c) shows the 
corresponding reduction of the line-width of the resonance peak. 

cycle, with Vp > 0. While the strain-dependence of the resonance freqnency, shown in panel 
(a), is weak, with Afo/fo ~ 0.02, a remarkably strong strain-dependence of the Q-factor 
is observed, with AQq/Qq ~ 0.25 over the same voltage range (panel (b)). The tnning 
cycle is repeated, and the dependence of fo and Qo are calculated for each compression and 
tension cycle, and collected in the histograms shown in the insets |Z7|. Fitting a Gaussian 
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distribution yields a mean frequency dependence of Afo = 6.3kHz/V and a Q-factor 
dependence of AQo = 0.09 /V. 

To calculate the induced strain as a function of the applied voltage, Vp, we consider 
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FIG. 3: (a,inset) optical image (left) and atomic force microscopy topography image (right) of 
device B. The drum is buckled, as is confirmed by the height profile shown in the main panel, 
(b) Qualitative behaviour of the resonance frequency in the pre-buckled and post-buckled regimes. 
Dotted: idealized (symmetric) resonator; solid: realistic (non-symmetric) device. (c,d) Measure¬ 
ments of the tuning of the resonance frequency and Q-factor when varying the compressive strain. 


a bimorph geometry |2H], with tsi and tp the thickness of the silicon and the piezo sheet. 
The respective Young’s moduli are Es, = 150 GPa and Ep = 62 GPa, and the thicknesses 
tsi = 100/im and tp = 127pm. With the piezoelectric coefficient hsi = —190 x 10“^^ 
and h = 0.534 a dimensionless number which represents the ratio’s of the Young’s moduli and 
the thicknesses, the strain in the M 0 S 2 is calculated as e — 53ih/tp X Vp — 8.0 X 10 V 
The calculated strain is plotted on a secondary x-axis in Fig.2. When the compressive 
strain exceeds a critical limit, the plate buckles. For a circular plate, the critical strain is 
calculated as = c’'cr/EMoS 2 = ES]- Here, EmoS 2 t^e real part of the Young’s 
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modulus, u = 0.25 is Poisson’s ratio, and t and r are the plate thickness and radius. K is 
a constant that depends on the boundary condition, with K = 1.22 for a clamped plate. 
Although in the present experiments the critical strain should occur at €„ = 4.7 x 10“^, 
which corresponds to Vp = —59 V, no buckling is observed 

We now turn our attention to device B, which is shown in detail in the inset of 
Fig.3(a). Clearly, a part of the drum is bulged: the bright color in the center indicates 
a buckle, which is the result of residual compressive strain which is introduced during 
fabrication. This device allows us to investigate the strain-dependent behaviour in the 
post-buckled regime. The topographic AFM image in the second inset, conhrms the 
presence of a buckle: a line-cut across the buckle (main hgure) reveals a height of several 
tens of nanometers, with multiple smaller corrugations and wrinkles superimposed. 

The presence of wrinkles make and analytical treatment of the resonance frequencies 
difficult, and instead, we present here only qualitatively the strain-dependent resonance 
frequency of a clamped-clamped plate. Figure 3(b) represents two cases: an idealized 
symmetric system (grey dots), and a system that more closely resembles the non-symmetric 
M 0 S 2 nanodrum (blue solid line). While in the pre-buckled regime the resonance frequency 
increases with the tensile strain, in the post-buckled regime it is expected to decrease 
with the tension. These situations apply to the devices considered: device A is in the 
pre-buckled regime and tunes according to the green arrow, and device B is post-buckled 
and is expected to tune along the blue arrow. 

Figure 3(c,d) show the measured tuning behaviour of drum B. Indeed, the frequency 
dependence is opposite to the one observed for drum A, as the resonance frequency 
decreases with the tensile strain. Compared to the flat device, the strain-dependence 
of the Q-factor is less pronounced. The low strain-dependence of the Q-factor could be 
explained by a relaxation of the compressive stress in the post-buckling regime, where 
elastic energy is converted from compression to bending. Interestingly, the frequency tuning 
curve shows a clear hysteresis: when the plate is compressed, the frequency vs. strain 
response follows a different path than when the strain is released. This cannot be due to 
hysteresis in the piezo stack, since in the measurements on device A in Fig. 2 the forward 
and backward tuning curves coincide. The observed effect is attributed to a change in the 
mechanical properties of the flake. Hysteresis in the post-buckled regime could indicate a 
conformational change, possibly of one of the wrinkles. Similar hysteretic effects could occur 
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at a smaller dimensional scale, and give rise to energy dissipation, causing the excessive 
damping of resonant motion. 

Another explanation for the strong strain-dependent damping could be the inevitable 
presence of (static) microscopic corrugations and wrinkles. Theoretical investigations have 
shown that microscopic geometric artefacts act as long-wavelength elastic scatterers EH. 
carrying away energy from the flexural modes. The wrinkles are not present in top-down 
silicon-based devices, which are inherently flat due to their fabrication process. Applying 
tensile strain to the 2-D resonator 'irons-out' the static wrinkles, which reduces the 
number of scatterers and results in a lower dissipation (i.e., a reduction of E 2 ), while the 
resonance frequency (Ei) is affected only weakly. In addition to the static wrinkles, the 2-D 
material resonators exhibit dynamic wrinkles due to the thermal fluctuations. Applying 
strain increases the spring constant (Ei), which reduces the mean squared amplitude of 
these fluctuation-induced dynamic wrinkles. Dynamic wrinkles are far less pronounced in 
top-down fabricated devices, which are typically thicker by one or two orders of magnitude 
and therefore have a much higher spring constant. This results in thermal fluctuations with 
relatively low amplitudes. 

Besides tuning the damping in mechanical resonators, there are other interesting 
applications for controlled strain tuning in 2-D materials. In these materials, which can 
be excessively strained due to the lack of defects [32], the mechanical strain changes the 
band structure. The qualitative changes in the electronic and optical properties can 

enable applications such as piezo-electric energy harvesters [511 EH] and pressure, motion, 
and mass sensors EZI- While biaxial strain can be adjusted by varying the temperature 
by deploying the thermal expansion mismatch [38[, the controlled application of strain 
described here, can be used to study the strain-dependent properties of 2-D materials in 
great detail. 

In conclusion, we studied experimentally the strain-dependence of the Q-factor in 
thin M 0 S 2 drum resonators in the pre- and post-buckling regime. The experiments indicate 
that, as in MEMS and NEMS resonators, the Q-factor increases with the applied tensile 
strain. However, in the M 0 S 2 resonators, the increase in Q is manifests as a reduction of 
the line-width, which indicates a decrease in dissipative part of the spring constant. This 
is in contrast to top-down fabricated MEMS resonators, where Q increases with strain due 
to an increase in the conservative part of the spring constant, which has only a small effect 


on the resonance line-width. This result sheds light on the very low Q-factors observed in 
recent experiments with 2-D mechanical resonators, and suggests that microscopic wrinkles 
and corrugations, which are ironed-out by applying tensile strain, could play a role in the 
observed low Q’s. In the post-buckled device, hysteresis is observed, and the Q-factor 
depends less on the strain. The experiment shows that strain engineering is a viable tool 
to reduce the damping in nanodrum resonators made from 2-dimensional materials. 
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Agreement no 318287, project LANDAUER. 


[1] J. S. Bunch, A. M. Van Der Zande, S. S. Verbridge, I. W. Frank, D. M. Tanenbaum, J. M. 
Parpia, H. G. Craighead, and P. L. McEuen, Science 315, 490 (2007). 

[2] C. Chen, S. Rosenblatt, K. I. Bolotin, W. Kalb, P. Kim, I. Kymissis, H. L. Stormer, T. F. 
Heinz, and J. Hone, Nature Nanotech 4, 861 (2009). 

[3] A. Eichler, J. Moser, J. Chaste, M. Zdrojek, I. Wilson-Rae, and A. Bachtold, Nature Nanotech 
6, 339 (2011). 

[4] R. A. Barton, B. flic, A. M. van der Zande, W. S. Whitney, P. L. McEuen, J. M. Parpia, and 
H. G. Craighead, Nano Lett 11, 1232 (2011). 

[5] J. Lee, Z. Wang, K. He, J. Shan, and P. X.-L. Feng, ACS Nano 7, 6086 (2013). 

[6] A. Castellanos-Gomez, R. van Leeuwen, M. Buscema, H. S. J. van der Zant, G. A. Steele, and 
W. J. Venstra, Adv Mater 25, 6719 (2013). 

[7] R. van Leeuwen, A. Castellanos-Gomez, G. A. Steele, H. S. J. van der Zant, and W. J. Venstra, 
Appl Phys Lett 105, 041911 (2014). 

[8] S. S. Verbridge, H. G. Craighead, and J. M. Parpia, Appl Phys Lett 92, 013112 (2008). 

[9] A. Croy, D. Midtvedt, A. Isacsson, and J. M. Kinaret, Phys Rev B 86, 235435 (2012). 

[10] A. Eriksson, D. Midtvedt, A. Croy, and A. Isacsson, Nanotechnology 24, 395702 (2013). 

[11] A. Laitinen, M. Oksanen, A. Fay, D. Cox, M. Tomi, P. Virtanen, and P. J. Hakonen, Nano 


9 



Lett 14, 3009 (2014). 

[12] D. Midtvedt, A. Croy, A. Isacsson, Z. Qi, and H. S. Park, Phys Rev Lett 112, 145503 (2014). 

[13] C. Edblom and A. Isacsson, Phys Rev B 90, 155425 (2014). 

[14] S. S. Verbridge, D. F. Shapiro, H. G. Craighead, and J. M. Parpia, Nano Lett 7, 1728 (2007). 

[15] Q. P. Unterreithmeier, T. Faust, and J. P. Kotthaus, Phys Rev Lett 105, 027205 (2010). 

[16] S. Y. Kim and H. S. Park, Nanotechnology 21, 105710 (2010). 

[17] Z. Y. Ning, T. W. Shi, M. Q. Fu, Y. Guo, X. L. Wei, S. Gao, and Q. Chen, Nano Lett 14, 
1221 (2014). 

[18] S. Schmid and C. Hierold, J Appl Phys 104, 093516 (2008). 

[19] S. Schmid, K. D. Jensen, K. H. Nielsen, and A. Boisen, Phys Rev B 84, 165307 (2011). 

[20] P.-L. Yu, T. Purdy, and C. Regal, Phys Rev Lett 108, 083603 (2012). 

[21] C. Chen, S. Lee, V. V. Deshpande, G.-H. Lee, M. Lekas, K. Shepard, and J. Hone, Nature 
Nanotech 8, 923 (2013). 

[22] J. S. Bunch, S. S. Verbridge, J. S. Alden, A. M. van der Zande, J. M. Parpia, H. G. Craighead, 
and P. L. McEuen, Nano Lett 8, 2458 (2008). 

[23] M. K. Zalalutdinov, J. T. Robinson, C. E. Junkermeier, J. C. Culbertson, T. L. Reinecke, 
R. Stine, P. E. Sheehan, B. H. Houston, and E. S. Snow, Nano Lett 12, 4212 (2012). 

[24] S. J. Cartamil, P. G. Steeneken, F. D. Tichelaar, E. Navarro-Moratalla, W. J. Venstra, van 
Leeuwen R, E. Coronado, H. S. J. van der Zant, G. A. Steele, and A. Castellano-Gomez, Nano 
Research, in press pp. http;//dx.doi.org/10.1007/sl2274-015-0789-8 (2015). 

[25] A. Castellanos-Gomez, M. Buscema, R. Molenaar, V. Singh, L. Janssen, H. S. J. van der Zant, 
and G. A. Steele, 2D Materials 1, 011002 (2014). 

[26] N. O. Azak, M. Y. Shagam, D. M. Karabacak, K. L. Ekinci, D. H. Kim, and D. Y. Jang, Appl 
Phys Lett 91, 093112 (2007). 

[27] The histograms represent 121 measurements cf. Fig. 2(a) for the frequeney tuning, and 56 
measurements ef. Fig.2(b) for the tuning of the Q-faetor. 

[28] S. Timoshenko, J. Opt. Soc. Am 11, 233 (1925). 

[29] W. C. Young and R. G. Budynas, Roark’s formulas for stress and strain, vol. 7 (McGraw-Hill 
New York, 2002). 

[30] No buekling was observed for voltages up to Vp = —150 V. 

[31] E. E. Helgee and A. Isacsson, Phys Rev B 90, 045416 (2014). 


10 



[32] H. H. Perez Garza, E. W. Kievit, G. F. Schneider, and U. Stanfer, Nano Lett 14, 4107 (2014). 

[33] K.-A. N. Duerloo and E. J. Reed, Nano Lett 13, 1681 (2013). 

[34] W. Wu, L. Wang, Y. Li, F. Zhang, L. Lin, S. Niu, D. Chenet, X. Zhang, Y. Hao, T. F. Heinz, 
et ah, Nature 514, 470 (2014). 

[35] B. Amorim, A. Cortijo, F. de Juan, A. Grnshin, F. Guinea, A. Gutierrez-Rubio, H. Ochoa, 
V. Parente, R. Roldan, P. San-Jose, et ah, arXiv preprint arXiv:1503.00747 (2015). 

[36] M. Lopez-Suarez, M. Pruneda, G. Abadal, and R. Rurali, Nanotechnology 25, 175401 (2014). 

[37] C. Yan, J. Wang, W. Kang, M. Cni, X. Wang, C. Y. Foo, K. J. Ghee, and P. S. Lee, Adv 
Mater 26, 2022 (2014). 

[38] W. Bao, F. Miao, Z. Chen, H. Zhang, W. Jang, C. Dames, and C. N. Lan, Nature Nanotechn 
4, 562 (2009). 


11 



